Calorimetric Evidence for Nodes in the Overdoped Ba(Feo.gCoo.i)2As 2 
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We present low-temperature specific heat of the electron-doped Ba(Feo.9Coo.i)2As2, which does 
not show any indication of an upturn down to 400 mK, the lowest measuring temperature. The 
lack of a Schottky-like feature at low temperatures or in magnetic fields up to 9 Tesla enables 
us to identify enhanced low-temperature quasiparticle excitations and to study anisotropy in the 
linear term of the specific heat. Our results can not be explained by a single or multiple isotropic 
superconducting gap, but are consistent with multi-gap superconductivity with nodes on at least 
one Fermi surface sheet. 

PACS numbers: 74.20.Rp, 74.25. Fy, 74.25.Bt 
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Since the discovery of the iron-based pnictide su- 
perconductors, substantial experimental and theoreti- 
cal work has been performed to understand the super- 
conducting (SC) mechanism and its manifestation via 
the gap structure. [l| Despite these efforts, however, 
the gap symmetry is still unclear. In RFeAsO (R: 
rare earth), the Till' phase, a full gap is expected 
from Andreev spectroscopy, [2| but a nodal gap is re- 
ported from optical measurements Q for a compound 
with similar stoichiometry. Conflicting interpretations 
of the SC order parameter also have been reported in 
Co-doped BaFe2As2 (Bal22 phase). Angle-resolved pho- 
toemission spectroscopy (ARPES) [ij and scanning tun- 
neling microscopy (STM) [B[ suggest an isotropic gap; 
whereas, Raman scattering, [(| nuclear magnetic reso- 
nance (NMR) Ml, m uon spin rotation (/nSR) Q , thermal 
conductivity [3, [TfJ, specific heat, [Tljj and penetration 
depth [l^, El measurements all indicate an anisotropic 
gap on one or more Fermi surface sheets. Strong doping 
dependence of superconducting properties and/or spe- 
cific features of the measurement techniques which probe 
preferentially a fraction of the total Fermi surface could 
be behind these seemingly contradicting results. [H], HH 

Low-temperature specific heat measurements, which 
probe the electronic density of states, have been in- 
strumental in determining the SC order parameter of 
various classes of superconductors, including high-T c 
cuprates. lB"18| In the newly discovered Fe- based pnic- 
tides, the interpretation of specific heat measurements at 
T <C T c has been controversial because of a Schottky- 
like upturn at the lowest temperatures, which has been 
attributed to some kind of impurities. 14, [I9I - I21I ] Here 



we report low-temperature specific heat measurements 
of Ba(Feo.gCoo.i)2As2, in which there is no Schottky- 
like feature down to 400 mK and up to 9 T, the lowest 
measured temperature and the highest applied magnetic 
field, respectively. The lack of Schottky contribution al- 
lows us to identify the low-T electronic excitations and 



to study the field dependence of the linear specific heat 
term. Our results are not compatible with an isotropic 
s-wave gap, but are consistent with gap zeroes on parts 
of a Fermi surface in a multiband model. 

Single crystalline Ba(Feo.gCoo.i)2As2 was grown using 
a closed Bridgman method with a tungsten crucible at 
1550 °C to contain volatile arsenic (As) [22}. Powder 
x-ray resonance diffraction shows that the lattice param- 
eters are a=3.9458 A and c=12.9476 A. These crystals 
have a plate-like tendency with the [001]-direction per- 
pendicular to the crystal plane. Specific heat was mea- 
sured down to 400 mK and up to 9 T in a Quantum 
Design PPMS (Physical Properties Measurement Sys- 
tem) with a 3 He option. The superconducting transi- 
tion temperature, defined from the midpoint of the spe- 
cific heat discontinuity (see Fig. 2a), is 19.0 K, which is 
lower than that of optimally doped Ba(Feo.92Coo.os)2As2 
(T c = 25 K) and indicates that the present compound is 
slightly overdoped. 

Figure 1 is a plot of the temperature-dependent spe- 
cific heat of Ba(Feo.9Coo.i)2As2 in zero field. A specific 
heat jump at 19.0 K signals the superconducting phase 
transition, confirmed by low-field magnetic susceptibility 
measurements. Unlike magnetic and structural transi- 
tions in Ba(Fei_ x Coa;)2As2, which are strongly affected 
by Co-doping, the phononic contribution to the specific 
heat is almost independent of doping. 23, [HI In the 



ensuing analysis, therefore, we use the phonon specific 
heat measured from undoped BaFe2As2 as the phonon 
contribution in the present material. Because there is 
no Schottky-like feature at the lowest temperatures of 
these measurements, the normal-state specific heat (C n ) 
can be written as C n /T — j n + C P h/T, where j n is 
the normal-state Sommerfeld coefficient and C p h is the 
phonon contribution. The solid curve in Fig. 1 is the nor- 
mal state C n /T with 7„=23.5 mJ/mol-K 2 , which satisfies 
entropy conservation at T c (see the inset of Fig. 1 and 
discussion below). In the superconducting state, there 
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FIG. 1: (color online) Specific heat divided by temperature 
of Ba(Feo.9Coo.i)2As2 in zero applied field (filled squares). 
The red solid line is the phonon contribution C p h/T deter- 
mined from the parent compound BaFe2As2, plus a constant 
7n=23.5 mj/mol K 2 that is required to balance entropy at 
T c . Inset: Electronic entropy of normal (red circles) and su- 
perconducting (black squares) states. A deviation of Sea from 
Sen appears near 19 K. 



is a noticeable T— linear specific heat contribution 70 of 
3.4 mJ/mol-K 2 , about 14 % of the normal-state j n . A siz- 
able 70 is often reported in superconductors with gap ze- 
roes on parts of the Fermi surface because impurities eas- 
ily produce a finite density of states even at zero energy. 
This finite value of 70 suggests the possibility of nodal 
superconductivity in Ba(Feo.gCoo.i)2As2. Even though 
a finite 70 could arise from non-superconducting parts of 
the crystal, the lack of a Schottky-like feature in C/T 
and lack of a second impurity phase in x-ray diffraction 
suggest that any contribution from impurities is negligi- 
bile in the ensuing analysis. Supporting this conclusion, 
the value of 70 is smaller than that reported earlier on 
a compound with similar Co concentration and an up- 
turn in the low— T specific heat (2^. We note that in 
multiband systems, a zero-energy impurity band can be 
effectively formed without nodes on any Fermi surface 
sheet, resulting in large value of 70 (26j . 

Electronic entropies of the normal (S en ) and supercon- 
ducting (S es ) phases of Ba(Feo.gCoo.i)2As2 are plotted in 
the inset of Fig. 1 and are used to estimate various SC 
properties. The superconducting condensation energy, 
U = J^ c (S en — S es )dT, is estimated to be 1.23 J/mol. 
The zero-temperature thermodynamic critical field B c (0) 
is 0.22 T, from the relation U = B 2 (0)/2^ - When 
combined with previous measurements of the penetra- 
tion depth (A = 325 nm) and coherence length (£ = 
27.6 nm), 



, [271 |28[ the upper critical field is esti- 
mated to be 37 T from the relation ^oH c2 = \/2kB c (0), 
k = A/£ = 11.8. These superconducting properties of 
the present sample with x = 0.10 are quite comparable 
to those reported previously on compounds with similar 
stoichiometry. 

Figure 2 displays the low-temperature electronic spe- 
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FIG. 2: (color online) a Temperature dependence of the 
electronic specific heat. The experimental data are com- 
pared with the two-band model described in the text, with 
an isotropic gap on the hole sheet, and a nodal gap on the 
electron sheet. Note that, when the nodal behavior is weak 
(r — 0.55, nodes close to each other), the numerical result 
using the scattering rates required to explain the residual lin- 
ear term in the superconducting state exhibits a downturn 
at low temperatures, not observed experimentally. Therefore 
the low-temperature data are most consistent with a relatively 
well-pronounced nodal behavior, r > 0.6. b Magnification of 
the electronic specific heat at low temperatures. 



cific heat of Ba(Feo.gCoo.i)2As2 after subtracting the 
phononic contribution. As seen in Fig. 2b, there is no 
hint of an upturn in C e i/T down to 400 mK and up 
to 9 T, irrespective of the direction of the applied field 
(see Fig. 3). Thus the closed Bridgman method used 
for preparing this crystal seems to avoid incorporation of 
impurity phases that are generated from a flux-growth 
method, enabling us to study the low-temperature spe- 
cific heat without the need to subtract poorly understood 
values to account for the field-dependent upturn in pre- 
viously studied samples. 

To investigate the gap symmetry, we begin by analyz- 
ing the temperature dependence of C e i/T in zero field, 
where C e i = C — C p h ■ We analyze the data using an ef- 
fective two band model with an isotropic gap on the hole 
Fermi surface and a gap on the electron Fermi surface of 
the form A e (k) = A e (l — r + rcos2^), where <j> is the 
angle measured from the [100] direction. This form is 
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FIG. 3: (color online) Low-temperature specific heat plotted 
for various magnetic fields applied parallel and perpendicular 
to the crystalline c-axis in panel a and b, respectively. Solid 
lines at zero field are guides to eyes that show a C oc T 2 
dependence. 



29, 30]. and is consistent 



suggested by several theories 
with the general Ai g symmetry of the order parameter in 
the Brillouin Zone. In this expression for the gap func- 
tion, r > 0.5 (r < 0.5) corresponds to the nodal (fully 
gapped) order parameter on the electron Fermi surface 
sheet. The gap magnitudes on the two Fermi surfaces 
were computed self-consistcntly assuming a purely in- 
terband pairing interaction requiring a sign change of 
the order parameter between the electron and the hole 
sheets 
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325. The residual linear, in T, term and the 
temperature variation of the electronic specific heat are 
controlled by the value of r and by disorder. In our model 
the impurities are characterized by their concentration, 
n impi an d the strength of the inter- and intra-band scat- 
tering, U\ and Uq respectively. The quasiparticle scat- 
tering rate T depends also on the density of states at 
the Fermi surface, Nf, which , for simplicity, we take to 
be the same on the two sheets. We use the transition 
temperature of the pure (impurity-free) sample, T c q, as 
a unit of energy, and hence introduce a dimensionless 



scattering rate r\ = T n /2nT c o, where F r 



,/TtNf. 



Together with the phase shift of the intraband scatter- 
ing, 5 = arctan(irNfUo), and the ratio of the inter- to 
intra-band potential SV = Ui/Uq it completely deter- 
mines the properties of the impurity ensemble. Since the 
inter- (intra-) band potential describes scattering with a 
large (small) momentum transfer, for usual charged im- 
purities we expect SV < 1, although SV « 1 is plausible 
for the direct substitution to the Fe site that strongly af- 
fects the d orbitals. Indeed, our best fits give a moderate 
ratio SV. We determine the quasiparticle self-energies 
using the self-consistent T— matrix approximation [33l — 
351 ] . and compute the density of states, the entropy, and 
the specific heat for comparison with the experimental 
results. 

We found that, irrespective of the strength of the 
scattering and the impurity concentration, the low- 



temperature data cannot be reliably fit using a model 
with a fully gapped electron Fermi surface, r < 0.5. Even 
though the density of states at the Fermi surface, and 
therefore the Sommerfcld coefficient hmj^ C e i(T)/T, 
may be finite in this case due to formation of an imp urity 
band in the presence of strong interband scattering [35| , 
more curvature is expected for the temperature depen- 
dence of C e i in that case than is observed experimentally. 
The nearly linear dependence in Fig. 2 can only be satis- 
factorily explained assuming a nodal gap structure, i.e., 
r > 0.5. Even in that case, the intraband scattering tends 
to lift the nodes [35[ , resulting in a pronounced downturn 
in C e i/T for the values of r close to the "critical" value 
that separates the gapless regime from that where the 
nodes are lifted by disorder and a full gap opens (see 
Fig. 2). As shown in Fig. 2, the best fits are for the 
well-defined nodes, r = 0.6. 

Measurements under applied magnetic field provide 
additional information on the gap structure. The elec- 
tronic specific heat divided by temperature C e i /T is dis- 
played against T for magnetic field perpendicular and 
parallel to the crystalline c-axis in Fig. 3a and 3b, respec- 
tively. With increasing field, the temperature variation of 
Cei/T at low T becomes weaker and almost disappears at 
9 T for H//c. Since the extrapolated T = intercept of 
C e i/T, 70, depends on the power law that we use for the 
extrapolation, in Fig. 4a we plot the normalized values of 
the C/T at 0.4 K (T/T c = 0.021) to study the magnetic 
field dependence of the density of states. For an s-wave 
SC gap, the electronic density of states is proportional 
to the density of vortices in the mixed state, leading to 
a linear field dependence of 70. In contrast, the presence 
of nodes in the SC gap changes the field-dependence to 
a sub-linear behavior due to a Doppler shift of the nodal 
quasiparticles. Supporting the conclusion of nodal super- 
conductivity from the low-T specific heat analysis, C e i/T 
which is proportional to the electronic density of states 
deviates from a linear- dependence on the field, but shows 
a sub-linear behavior. Comparison between the experi- 
mental data for field along the c-axis and the theory with 
SC gap nodes is shown in Fig. 4a, where we use the val- 
ues of r along with the impurity parameters from Fig. 2 
to compute the field dependence of the specific heat un- 
der the applied field using the modified Brandt-Pesch- 
Tewordt (BPT) approximation 36[. This approximation 
is designed for intermediate and high fields and is known 
to give only a qualitatively correct picture for nodal gaps 
at low H . As seen in Fig. 4a, the calculations reproduce 
the sublinear field dependence, consistent with the exis- 
tence of nodes. At the same time the magnitude of the 
field enhancement is much greater in theory than in ex- 
periment. This is due to a significant overestimate of the 
contribution of the fully gapped hole Fermi surface to the 
density of states within the BPT approximation at low 
fields 01[32|. 



The field dependence of the specific heat coefficient 
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FIG. 4: (color online) a Field dependence of the linear term 
in the heat capacity at 0.4 K (T/T c = 0.021) divided by the 
normal state Sommerfeld coefficient 7„. Theoretical calcu- 
lations using the same parameters that give the best fit to 
the temperature dependence of C e i/T in zero field reproduce 
the overall field dependence but yield a higher magnitude of 
the effect: the dashed line is calculated from the parameters 
(r = 0.6,5 = 80°, r] = 0.1, SV = 0.6), while the dash-dotted 
line is from (r = 0.6,5 = 80°, t) = 0.05,5V = 0.8). As ex- 
plained in text, the larger absolute value of the calculated 
y(H) is due to an overestimate of the contribution of the hole 
band within the modified Brandt-Pesch-Tewordt approxima- 
tion. For reference, we show least-squares fits of the data for 
H || c and H J_ c to a power law form a + bH* (dotted lines), 
where n = 0.7 reasonably explains the data, b Anisotropic 
field dependence of the linear term in the heat capacity at 
0.4 K. The ratio of the Sommerfeld coefficient for field par- 
allel and perpendicular to the crystalline c-axis, S^c/S^ab, is 
plotted against magnetic field, where S*y = *y(H) — 7(0T) 



(-f(H)) is fit by a power-law form a + bH 7 for both field 
directions with b = 0.64 and 0.33 for H \\ c and fflc, 
respectively (dotted lines in Fig. 4a). Anisotropy of the 
field-induced quasiparticle density of states, Sjc/Sjab, is 
plotted in Fig. 4b, where £7 = j(H) - j(0T). This 
anisotropy ratio is approximately two over the measured 
field range and reflects the superconducting anisotropy 
through 5-y c /5j ab = {H^/H^f 0. When e = 0.7 is 
used, the uppercritical field anisotropy ratio is approxi- 
mately 2.7. This value is compatible with that estimated 
from the initial slope (rfJLWrfT) at T c by the Werthamer- 
Helfand-Hohenburg model 39,40]. Kano et al., however, 



reported that the H C 2 anisotropy ratio at T = K in a 
compound with similar stoichiometry is close to 1, where 
the difference has been interpreted as due to multi-band 
effects in the pnictides 4l|. Indeed, there is substantial 
evidence for multiple electronic bands in the pnictides 
and their consequences for superconductivity. Studies of 
the Fermi surface find a distinct asymmetry between hole 
and electron bands: a circular hole Fermi surface is ob- 
served around the T point, but an anisotropic electron 
Fermi surface exists around the M point [3l|, |42T - |44T | . Re- 
cently, penetration-depth and thermal conductivity mea- 
surements have suggested that line nodes may exist on an 
electron sheet near the Brillouin zone center and that hole 



bands centered around the T point are fully gapped [45- 

Si. 

In conclusion, the specific heat of Ba(Feo.gCoo.i)2As2 
prepared by the closed Bridgman method does not show 
any hint of an upturn at low temperatures, a materi- 
als problem that has plagued crystals prepared by other 
methods. The lack of a Schottky-like feature enables 
us to obtain the intrinsic low-temperature specific heat 
whose temperature dependence is consistent with an or- 
der parameter that produces an isotropic SC gap on one 
of the Fermi surface sheets and a gap with line nodes 
on the second, presumably electron, sheet. Evidence for 
nodal superconductivity is supported by the magnetic 
field dependence of the Sommerfeld coefficient 7 and a 
finite residual density of states in the superconducting 
state. Additionally, these experiments suggest that a 
complete description of superconductivity requires the 
involvement of multiple bands. 

This work was supported by the National Research 
Foundation (NRF) grant (2010-0016560) funded by Ko- 
rea government (MEST). Work at Los Alamos was per- 
formed under the auspices of the U. S. Department of 
Energy/Office of Science and supported in part by the 
Los Alamos LDRD program. ZY and TP acknowl- 
edges support by the promotion program for new fac- 
ulty, Sungkyunkwan University (2009). I.V. is supported 
in part by US DOE via Grant DE-FG02-08ER46492. 
ABV's support in part comes from US NSF grant DMR- 
0954342. Y.S.K. is supported by the Basic Science Re- 
search Program (2010-0007487) and Nuclear R&D Pro- 
grams (2006-2002165 and 2009-0078025). 



[1] K. Ishida, Y. Nakai, and H. Hosono, J. Phys. Soc. Jpn. 

78, 062001 (2009). 
[2] T. Y. Chen, Z. Tesanovic, T. H. Liu, X. H. Chen, and C. 

L. Chien, Nature 453, 1224 (2008). 
[3] A. Dubroka et al., Phys. Rev. Lett. 101, 097011 (2008). 
[4] K. Terashima et al., Proc. Natl. Acad. Sci. USA 106 7330 

(2009). 

[5] Y. Yin et al., Phy. Rev. Lett. 101, 097002 (2009). 

[6] B. Mushuler et al., Phys. Rev. B 80, 180510(R) (2009). 



5 



[30] 



F. Ning et al., J. Phys. Soc. Jpn. 77, 103705 (2008). 
T. J. Williams et al., Phys. Rev. B 80, 094501 (2009). [31 
M. A. Tanatar et al, Phys. Rev. Lett. 104, 067002 
(2010). [32 
J. K. Dong et al., Phys. Rev. B 81, 094520 (2010). 
Gang Mu et al., Chin. Phys. Lett. 27, 037402 (2010). [33 
R. T. Gordon et al., Phys. Rev. Lett. 102, 127004 (2009). 
K. Hashimoto et al., Phys. Rev. B 81, 220501(R) (2010). [34 
T. M. McQueen et al., Phys. Rev. B 79, 01452 2(2009). 
J. S. Kim et al., unpublished (|arXiv:1002.3355H L [35 
K. A. Moler, D. J. Baar, J. S. Urbach, R. Liang, W. 
N. Hardy, and A. Kapitulnik, Phys. Rev. Lett. 73, 2744 [36 
(1994). 

T. Park et al., Phys. Rev. Lett. 90, 177001 (2003). [37 
T. Sakakibaraet al., J. Phys. Soc. Jpn. 76, 051004 (2007). [38 
K. Gofryk et al., New. J. Phys. 12, 023006 (2010). 
Gang Mu et al., Phys. Rev. B 79, 174501 (2009). [39 
J. S. Kim, E. G. Kim and G. R. Stewart, J. Phys.: Con- 
dens. Matter 21, 252201 (2009). [40 
Y. J. Song et al., Appl. Phys. Lett. 96, 212508 (2010). [41 
N. Ni et al., Phys. Rev. B 78, 214515(R) (2008). [42 
J.-H. Chu, J. Analytis, C. Kucharczyk, and I. R. Fisher, 
Phys. Rev. B 79, 014506 (2009). [43 
K. Gofryk et al., Phys. Rev. B 81, 184518 (2010). [44 
Y. Bang, Phys. Rev. Lett. 104, 217001 (2010). [45 
L. Luan et al., Phys. Rev. B 81, 100501(R) (2010). [46 
C. Bernhard et al., New J. Phys. 11, 055050 (2009). [47 
A. V. Chubukov and I. Eremin, Phys. Rev. B 82, [48 
060504(R) (2010). [49 
T. A. Maier, S. Graser, D. J. Scalapino, P. J. Hirschfeld, 



Phys. Rev. B 79, 224510 (2009). 

I. I. Mazin, D. J. Singh, M. D. Johannes, and M. H. Du, 
Phys. Rev. Lett. 101, 057003 (2008). 
A. V. Chubukov, D. V. Efremov, and I. Eremin, Phys. 
Rev. B 78, 134512 (2008). 

P. J. Hirschfeld, P. Wolfle, and D. Einzel, Phys. Rev. B 
37, 83 (1988). 

A. V. Balatsky, I. Vekhter, and J.-X. Zhu, Rev. Mod. 
Phys. 78, 373 (2006). 

V. Mishra, A. Vorontsov, P. J. Hirschfeld, and I. Vekhter, 
Phys. Rev. B 80, 224525 (2009). 

A. B. Vorontsov and I. Vekhter, Phys. Rev. B 75, 224501 
(2007). 

W. Pesch, Z. Phys. B 21, 263 (1975). 

M. Ichioka, A. Hasegawa, and K. Machida, Phys. Rev. B 

59, 184 (1999). 

N. R. Werthamer, E. Helfand, and P. C. Hohenberg, 
Phys. Rev. 147, 295 (1966). 

Y. Mamoto et al., Appl. Phys. Lett. 94, 062511 (2009). 
M. Kano et al., J. Phys. Soc. Jpn. 78, 084719 (2009). 
K. Seo, B. A. Bernevig, and J. Hu, Phys. Rev. Lett. 101, 
206404 (2008). 

Rev. B 80, 024515 (2009). 
Rev. Lett. 101, 177005 (2008). 
Phys. Rev. Lett. 102, 147001 (2009). 
Phys. Rev. Lett. 103, 127003 (2009). 
Phys. Rev. B 80, 220509(R) (2009). 



M. Yi et al., Phys. 
C. Liu et al. Phys. 
J. D. Fletcher et al. 
C. W. Hicks et al., 
M. Yamashita et al.. 



K. Hashimoto et al., arXiv:0907.4399 (unpublished) 
J.-Ph. Reid et al., Phys. Rev. B 82, 064501 (2010). 



